Weak gravitational lensing by fourth order gravity black holes 
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We discuss weak lensing characteristics for black holes in a fourth order / (R) gravity theory, 
characterized by a gravitational strength parameter a and a distance scale r c . Above r c gravity is 
strengthened and as a consequence weak lensing features are modified compared to the Schwarzschild 
case. We find a critical impact parameter (depending upon r c ) for which the behavior of the 
deflection angle changes. Using the Virbhadra-Ellis lens equation we improve the computation of 
the image positions, Einstein ring radii, magnification factors and the magnification ratio. We 
demonstrate that the magnification ratio as function of image separation has a different power-law 
dependence for each parameter a. As these are the lensing quantities most conveniently determined 
by direct measurements, future lensing surveys will be able to constrain the parameter a based on 
this prediction. 
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I. INTRODUCTION 
t-H ■ 

■ The recent advent of the so-called "Precision Cosmology" along with galactic observations indicate that General 
| Relativity (GR) with standard matter sources disagrees with an increasing number of observational data, e. g. those 

00 . coming from IA-type Supcrnovae, used as standard candles, large scale structure ranging from galaxies up to super- 
clusters [IHHj and galactic rotation curves. In addition, from a theoretical point of view, being not renormalizable, 
GR fails to be quantized in any standard way (see Q). Therefore from ultra-violet up to infrared scales GR is not 
^~-) ■ and cannot be the definitive theory of Gravitation despite the fact that it successfully addresses a wide range of 
£SJ . phenomena and the Newtonian weak field limit is correctly recovered. 

In order to interpret the recent observational data in the framework of GR, the introduction of unknown dark matter 
(DM) (to address dynamical phenomena as the formation of self-gravitating astrophysical structures) and dark energy 
• £h ■ (DE) (to address the problem of cosmic acceleration) seems to be necessary: however, the price of preserving the 
' simplicity of the Hilbert-Einstein Lagrangian has been the introduction of rather odd-behaving physical entities that, 

■ up to now, have not been revealed by any experiment at fundamental scales. This situation has led to several attempts 
devoted either to recover the validity of GR at any scale, or to construct alternative gravity theories that suitably 



generalize the Einstcinian one. The philosophy of these two schemes is that in the former case one has to modify the 
matter sector introducing DM and DE, in the latter approach the dynamics of the geometry (i.e. the left-hand-side 
of the Einstein equations) is modified but with the constraint to recover GR at local scales. 

Higher-order theories of gravity (both in metric 0-0 and Palatini flfjl - [l2| formulations) represent an interesting 
approach able to fruitfully cope with both dark matter and dark energy problems. A further approach is based on 
scalar - tensor theories of gravity but it can be shown that higher-order theories and scalar tensor ones can be related 
by conformal transformations (see, e.g., @, [H[ and references therein). 

It has been widely demonstrated that such theories can agree with the cosmological observations of the Hubble flow 
[3 [TH and the large scale structure evolution [jq . In addition, in the weak field limit, the gravitational potential 
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turns out to be modified [17H20l | in such a way that interesting consequences to galactic dynamics may be achieved 
without violating, at the same time, the constraints on the parametrized post-Newtonian (PPN) parameters coming 
from Solar System tests [Il| . 

If alternative theories of gravitation are able to explain both cosmological and local observations without the intro- 
duction of exotic energy-momentum sources, then one might ask how would the differences between these alternative 
theories be compared to GR with the unusual sources? It is proposed that gravitational lcnsing might be able to act a 
means for determining which theory governs the gravitational interaction, through measurements of image separations 
and the brightnesses of those multiple images. 

It is well known that the deflection of light observed during the Solar eclipse of 1919 was one of the first experimental 
confirmations of Einsteinian GR. Gravitational lensing, i.e. the deflection of light rays crossing the gravitational field 
of a compact object referred to as the lens, has become one of the most astonishing successes of GR and it represents 
nowadays a powerful tool capable of putting constraints on the dynamics of gravitational structures at different scales, 
from stars to galaxies and clusters of galaxies, and from the large scale structure to cosmological parameters (22l. l23j . 
If we modify the Lagrangian of the gravitational field, it is obvious that gravitational lensing should be affected. It 
is therefore mandatory to investigate how gravitational lensing behaves in the framework of alternative theories of 
gravity to develop a further check on the existence of DM and DE. 

In particular, one has to verify that the phenomenology of standard gravitational lensing is recovered in the limit 
as the modified theory of gravity reduces to GR, since several observations point to the validity of GR. However, it 
is worth stressing that the presence of DM has to be invoked in such cases, in particular for large-scale structure (see 
e.g. the case of Bullet Cluster [HI)- On the other hand, it is worth exploring whether deviations from the classical 
results for the main lensing quantities could be detected and act as clear signatures for modified theories of gravity. 

As a first step towards such an ambitious task we consider here power -law fourth order theories, i.e. we replace the 
Ricci scalar R in the gravity Lagrangian with the function f(R) oc R n (with n a positive integer ) an d investigate how 
this affects the gravitational lensing in the case of black hole solutions to the theory. In Ref. [25[ the gravitational 
lensing in / (R) theories with a Yukawa-type correction in the potential was investigated and the conclusion reached 
was that weak lensing could not discern between these theories and general relativity. As an alternative, we will 
consider the weak lensing phenomenology for black hole space-time coming from i?"-gravity in which gravitational 
lensing can occur. 

The paper is organized as follows. The next section provides a short summary of f(R) theories, focusing on the 
special R n case. We will prove in Section Mil that the predictions of weak gravitational lensing are different in 
the fourth order theory and in general relativity. For this we determine the image locations, Einstein ring radii, 
magnification factors and the flux ratio, for various model parameters. We then prove that the flux ratio as function 
of image separation has a different power-law dependence for each model parameter. We summarize our findings in 
the Concluding Remarks. 



II. GENERALITIES OF /(i?)-GRAVITY 

Let us consider now a particular class of higher order theories of gravity, the /(i?)-gravity theories, that are the 
most natural extension of GR where the Hilbcrt-Einstcin action is modified with a general function of the scalar 
curvature 0,0- A general action is: 

A = J d i x^[f(R) + C m ] (1) 

where f(R) is a generic function of the Ricci curvature scalar R, differentiable at least up to the second order, g is 
the determinant of the metric and C m is the standard matter Lagrangian 1 Varying the action with respect to the 
metric components g^ v , one obtains the generalized field equations that can be recast as 0, H^] : 

i f i "I tH 

GW = jy^-A 2^ lf(R) - Rf'(R)} + f'Wwu - g,,of(R) \ + (2) 



It is possible to take into account also the Palatini approach in which the metric g and the connection T are considered independent fields 
(see for example [26p ) . Here we consider the Levi-Civita connection and use the metric approach. See 0, [27| for a detailed comparison 
between the two pictures. 
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where = i? M „ — — g^ v and T^ v are the Einstein tensor and the standard matter stress - energy tensor, respec- 
tively. The prime denotes derivative with respect to R. The two terms f'{R).„ v and Df'(R) imply fourth order 
derivatives of the metric g^ v so that these models are also referred to as fourth order gravity. Starting from Eq. §Z§ 
and adopting the Robertson - Walker metric, it is possible to show that the Friedmann equations may still be written 
in the usual form provided that an effective curvature fluid (hence the name of curvature quintessence) is added to 
the matter term with energy density and pressure depending on the choice of f(R) @. 

The simplest choice for f(R) is a power-law like / (R) oc R n with n the slope of the Lagrangian (clearly, with n = 1, 
we recover the Einstein theory). In order to compare these power-law models with GR, we can search for spherically 
symmetric solutions of Eqs. ([2J and then compute the Newtonian limit. One has to consider the gravitational held 
generated by a pointlike source and then solve the vacuum case. Under the hypotheses of a weak gravitational field 
and slow motion sources. In general, we can write the space-time metric as : 



ds 2 = A(r)dt 2 - B{r)dr 2 - r 2 dQ 2 (3) 

where dfl 2 — d9 2 + sin 2 Odip 2 is the line element on the unit sphere. It is worth noticing that the Birkhoff theorem 
does not hold for any f(R) but only for restricted classes of these models as discussed in details in (28j . 

To obtain the value of the two unknown functions A(r) and B(r), we take the 00- vacuum component and the trace 
of the field equations (J5J in the absence of matter : 



to get a single equation : 



3af(R)+Rf(R)-2f(R)=0 



l„ m J'(R);00 



f'(R) [3—-R + -/(#)- 3^^ = 0. (4) 

Eq. Q is completely general and holds for any function f(R). Note that even if the metric is stationary so that 
dt9^v = 0j the term /'(i?) ; oo is n °t vanishing because of the existence of non-zero Christoffcl symbols entering the 
covariant derivatives. With the above power-law choice for the /(i?)-Lagrangian, Eq. (Q| reduces to: 

2n-l . . . n - 1 dA (r) din R(r) 

while the trace equation reads : 



□i?"" 1 ^) = —lR n { r ) . (6) 
3n 

Note that as soon as n = 1, Eq. ([6]) reduces to R = 0, which, when inserted into Eq. ([5]), gives i?oo = and then 
the standard Schwarzschild solution is recovered. In general, expressing i?oo and R in terms of the metric ((31), Eqs. 
([5]) and (|6|) become a system of two nonlinear coupled differential equations for the two functions A(r) and B(r). A 
physically motivated hypothesis is to search for (Schwarzschild-like) solutions of the form 



A(r) 



2$(r) 



B(r) 



(7) 



where (f>(r) is the gravitational potential at the distance r from a pointlike mass m. 

With these hypotheses, the vacuum field equations reduce to a system of two differential equations in the only 
unknown function $(r). To be more precise, we can solve Eq. ([5]) or © to determine and then use the other 
relation as a constraint to obtain solutions of physical interest. 

A general exact black hole solution is 



<P{r;a,r c ) = — — 
Ir 



, r 
1+ I — 

Tc. 



(8) 
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FIG. 1: The ratio of the power-law f(R) gravitational potential and the Newtonian gravitational potential as a function of 
a £ [0, 1) and of the logarithm of the distance r > r c , normalized to rs (the Schwarzschild radius of the lens). We have chosen 
a typical bulge mass of 10 10 M Q es 2 x 10 40 kg and r c « r bu i ge « 10 20 m. 

where the gravitational potential differs from the standard Newtonian one due to the second term on the right hand 
side p9l |. containing the parameter a and a characteristic distance r c . It is straightforward to see that, for a = 0, 
the Newtonian potential is recovered and the metric reduces to the Schwarzschild one (GR-case). We exclude the 
cases er > 1 (as the correction to the Newtonian potential asymptotically diverges) and a = 1 (as the correction is a 
constant, obstructing asymptotic flatness). 

To check the validity of solution ©, we insert the expression for $(r) into Eqs. ([5]) and ©, which both reduce to 
algebraic equations and are solved to obtain : 



(n-l)(<r-3) [-a{l + a)V x r) a -*} 
where rj = r/r c , Vy — Gm/c 2 r c and 



aViVo 

VlT] 



v%v = o, 



(9) 



Vo = 3(cr - 3Yn A - {h£ - 31ct + 48)n z - (3^ - I6cr + 17)n - {a z - 4cr - 5) 



(10) 



Pi = 3(cr - 3) 2 (1 - a)n 3 + {a - 3) 2 (5cr - 7)n 2 - (3cr 3 - 17a 2 + Ma ~ 36)n + {a 2 -3a- A)a 



(11) 



Eq. @ is identically satisfied for particular values of n and a. It should be noted that when deriving Eq. ([5]) from 
Eq. we assumed R =/= 0, which eliminates the case n = 1. Secondly, a = 3 may also be rejected since as noted 
above, all a > 1 lead to divergences. The third factor cannot vanish as it would imply either a — or a = —I for 
which n is less than unity. 

We can then look for further solutions of Eq. ((9J solving 



Pi(n, a) V + aVMn, a) ~ ^(n, ct)?? = 



(12) 



We neglected the second term of Eq. (fT2|) since for large scale-lengths r c the parameter Vi <C 1 and both n and cr are 
of order unity. Eq. (TT2]) is an algebraic equation for a as function of n with the following three solutions : 
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3n-4 



12n 2 -ln-l - s/pij 
q(n) 

Yin 2 - 7n - 1 + y/p(i 
q(n) 



(13) 



with : 



p{n) = 36?i 4 + 12n 3 - 83n 2 + 50n + 1 , 

and 



q(n) = 6n 2 - An + 2 . 

It is easy check that, for n = 1 the second solution (fT3| gives <r = 0, i.e. the second approximate solution reduces to 
the Newtonian one as expected. As a final check, we have to insert back into the vacuum field Eqs. (fj| and © the 
modified gravitational potential (JS)) with: 




FIG. 2: The dependence of correction parameter a in the potential upon the exponent n. 

In the case n = 1 (implying a = 0) the potential reduces to the Newtonian $j\r, as expected. The potential also 
reduces to the Newtonian value at r — r c . For smaller values of r gravity is weakened compared to the Newtonian 
values, while for r > r c gravity is strengthened. 

While the power a of the correction term is a universal quantity (since it depends on the exponent n entering 
the gravity Lagrangian), the scaledength r c is related to the integration constants that have to be set to solve the 
fourth order differential equations of the theory. The radius r c can be considered as a further gravitational radius 
complementing the Schwarzschild radius, originating in the fact that we consider a fourth order theory (compared to 
the second order GR) and as such it introduces two further degrees of freedom of the gravitational field. We expect r c 
to be related to the peculiarities of each gravitational system. Therefore it can take different values depending upon 
the system's mass and typical length scale. 

The fact that gravity is strengthened above r c is illustrated in Fig. 1 by plotting the ratio of the potentials 
$ (r; cr, r c ) /$tv (?*) as function of a and r for the ranges < a < 1 and r > r c with a suitably chosen value of r c . In 
the case of a typical spiral galaxy we identify r c with the bulge radius rb u i ge . If we choose, for example, mass of the 
bulge as 10 10 M Q sa 2 x 10 40 kg, we find r c r bulge w 10 20 m. Therefore gravity is strengthened outside r hvlnp as 
compared to the Newtonian case, providing an alternative to dark matter as a source for a flat rotation curve [301 ] . 
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III. WEAK LENSING IN FOURTH ORDER GRAVITY 
A. The lensing geometry 



/ 




FIG. 3: The lensing geometry (see text for discussion of the symbols). Positive angli 
directed arcs. 



are represented with counterclockwise 



The lensing geometry is shown on Fig. |3l The optical axis OLN is defined by the observer position O, the lens 
position L , and intersects the source plane at N. In the source plane S represents the location of the source and 
Ii : 2 the locations of the two images. We use the notations Di s and D\ for the projection of the lens-source and 
observer-lens distances onto the optical axis respectively [3l[. The observer-source distance is D s = Di + Di s - The 
source is located at an angle /3 from the optical axis, chosen to be "above" the optical axis. Images are located at 
angles 9 1^2 with respect to the optical axis and they can be either positive (for the image above the optical axis) or 
negative (for the image below the optical axis). For either of the images we denote s =sgn (9, such that \9\ = s9. 
We follow the convention that the deflection angle is 6 > whenever the light is bent towards the optical axis, cf. 
Rcf. [33| . Similarly as in Ref. [3l| . we characterize the mass by the dimensionless parameter e = Gm/c 2 L, with 
L = D s D l /D ls . 

In Ref. [29| the leading order lens equation 



\O\-80- 



El 



-8 = 



(15) 



was employed for the discussion of the weak lensing. We increase the accuracy of the weak lensing analysis by applying 
here the Virbhadra-Ellis lens equation 



tan \9\ - tan (s/?) - — [tan \9\ + tan (S - \9\)} = 



(16) 



We compare the two approaches by computing the difference between the Einstein angles (defined as 9\ = 92, occurring 
when /3 = 0) obtained from Eqs. (| 1 5[) and (|16|) . respectively, as a function of the fourth order black hole parameter a 
and mass parameter e. The results are shown in Fig 2] 



B. The deflection angle 



In Ref. weak lensing by point-like / (R) black holes characterized by the potential ((8]) was investigated. A 
deflection angle 



6 = 



2Gm 
c 2 b 



v/tF(1 - <t)T(1 - Q-/2 ) ( b 
2r(3/2-cr/2) 



(17) 
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FIG. 4: The difference between the Einstein angles Be, obtained from Eq. (I16[) and 9eo, obtained from Eq. (|15[) . as function 
of the parameters a and e. We set the distances Di = 1 Mpc and D; s = 2 Mpc. 



was derived, with b = 8Di the impact parameter (defined as the distance of the lensing object to the straight line 
trajectory, which would occur in the absence of the lensing object). For a = the deflection angle reproduces the 
Schwarzschild value Sg = 4Gm/c 2 b, while in the limit a — >• 1 the deflection angle is one half of 6$. 

In order to investigate the behavior of 5 in between the limiting values, we represent the deflection angle as function 
of (7 in Fig. [5] Three conclusions stand out: i) by increasing the impact parameter at any fixed value of a, the deflection 
angle always decreases, as in the Schwarzschild case; ii) there is a critical value of the ratio b/r c at (b/r c ) crit = 2, below 
which the deflection angle monotonically decreases with increasing a, and above which there is a single maximum at 
some cr,5 max ; iii) the parameter value <r<5 max increases with the value of b/r c . The rate of decrease of 6 with increasing 
b is lessened as compared to the Schwarzschild case for small impact parameters of order r c . 




FIG. 5: The plots show the 5 (a) dependence for different values of b/r c , in units 2Gm/c r c = 1. The respective values of b/r c 
are from top to bottom 1 (dashed), 2 (solid), 4 (diamond) on panel-a and 10 (cross), 100 (circle), 1000 (dotted) on panel-b. 
The critical behavior appears at b/r c = 2. 



C. Image positions and magnifications 

In Ref. [2^ | the image positions, the radius of the Einstein ring, image magnifications and the Paczynski curve in 
microlensing experiments were also estimated. In this follow-up paper we improve on the accuracy of the weak lensing 
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FIG. 6: The image positions 8 as function of e and /3 for a = 0.25 (left) and a = 0.75 (right), for the distances Di = 1 Mpc 
and Die = 2 Mpc. The angle (3 is varied up to 0.0015 rad, similarly as on Fig 4b of Ref. [3j|. With decreasing (3, the image 
separations shrink accordingly. At (3 = the angle 8 represents the angular radius of the Einstein ring. As we expect the (3 = 
sections of the surfaces are symmetric with respect to the plane 8 = 0. 



characteristics, by employing a more accurate lens equation with the focus on the behavior of dimensionless observable 
quantities which can be derived from the image positions and magnification factors. Our task is to determine the 
measurable differences between the predictions of the fourth order theory and those of GR. 

A weak lens equation was derived in Ref. [3ll | exclusively by trigonometric considerations and was applied to the 
computations of the image positions, magnifications and flux ratios to second order accuracy (both in the mass-related 
and tidal charge related small parameters) for brane-world black holes. It has been also shown, how the Virbhadra- 
Ellis lens equation follows as an approximation (agreement is reached in the first order of the perturbations). For the 
purposes of the present paper, not being concerned with second order effects, we will employ the Virbhadra-Ellis lens 



equation |32J- 34[, together with the deflection angle (jTTJ) derived in Ref. [29|. This generalizes the approach of Ref. 
[29j, where the leading order lens equation [Eq. (12) there] was employed for the discussion of weak lensing effects. 

The numerical solution of the system of equations (fT7|) and (fl6]) gives the positions of the images as function of e 
and (3, represented on Fig. [6l for a = 0.25 and 0.75, respectively. In both cases decreasing (3 decreases the image 
separations. As expected, the (3 = sections give symmetric curves with respect to the planes 8 = 0. This is because at 
(3 = the angle 8 = 8e represents the angular radius of the Einstein ring. For small impact parameter (implying small 
(3) the image separations obey (6>i — 02) a=o 25 < (#i — ^2) CT=0 75 , whereas for large values of (3 the image separations 
behave as (8\ — #2)^=025 > i® 1 ~ ^^)a=0 75- This indicates that our analysis based on the Virbhadra-Ellis lens 
equation is more accurate than the first post-Newtonian order calculation performed in Section VIII. in Ref. [38j |. 
which states that f(R) gravity is indistinguishable from general relativity and is consistent with the observational 
value of the post-Newtonian parameter 7 = 1 + (2.1 ± 2.3) x 10~ 5 . 



The magnification of the images are defined as 

d8i,2 812 



Mi. 2 



d(3 (3 



(18) 



For Schwarzschild lensing /ii — /i2 = 1 always holds. Fig [7] shows the image separations and magnifications as 
functions of (3/8e for a = 0.25 and 0.75. The upper and lower solid curves represent the primary and secondary 
image magnification factors, respectively; their ratio is the dashed curve; and the dotted curve is the image separation 
normalized with respect to the Einstein angle. The strongest effect appears on the ratio of magnifications, which 
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FIG. 7: The image separations and magnifications as functions of /3/0e for a — 0.25 (left) and a = 0.75 (right). We fixed 
e — 3.375 x 10 -9 , while the distances Di = 1 Mpc and Di s — 2 Mpc were chosen for the plots. The upper and lower solid 
curves represent the primary and secondary image magnification factors, respectively; their ratio is the dashed curve; and the 
dotted curve is the normalized image separation. 



TABLE I: The different slopes k of the curves indicate power law behaviours with different exponents for the various a. 



(T 





0.25 


0.5 


0.75 


0.9999 


K 


6.8 


7.3 


7.55 


7.7 


8.1 



increases more rapidly with (3 as a increases. The reason for this is that the primary image is magnified more for 
larger values of a. 



D. Power-law behavior 



In a lensing observation the most straightforward measurements are 1) the angular separation between the two 
images and 2) the ratio of the magnification factors. The first measurement does not require information on the 
location of the lens position, needed to define the individual image positions. The second measurement does not 
require an absolute measure of image brightness, since we arc taking a ratio between the two. 

Therefore we plot the ratio of the magnification factor for the primary and secondary images as function of the image 
separation (normalized to the Einstein angle), on the log-log scale, in Fig [5] The black (lowest) curve characterizes 
the lensing by a Schwarzschild black hole (a = 0), the colored curves, from bottom to top correspond to the fourth 
order gravity black hole lensing for the parameter values a = 0.25, 0.5, 0.75 and 0.9999. 

Similar to Ref. [3l| . for image separations greater than about 2.5 times the Einstein angle, the ratio of the 
magnification factors for each value of a obeys a power-law relationship 



V 0e 

The different slopes n of the curves indicate power- law behaviors with different exponents, which are presented in 
Table H 

Given a large enough number of measurements of image separations and image brightnesses, these power-law 
relations provide an observational signature that can distinguish which among the fourth order f(R) theories with 
different a (or n) describes the black hole space-time. 
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<s=0.99»9 
n=0.75 



a=0.2S 



6,-6; 



0.4 



log 



Q E 



FIG. 8: The ratio of the magnification factor of the primary and secondary images as function of the image separation normalized 
to the Einstein angle, on the log-log scale for a = 0., 0.25, 0.5, 0.75, 0.9999. We fixed e = 3.375 x 10 _9 .We set the parameters 
Di — 1 Mpc, Di s = 2 Mpc. The different slopes of the curves indicate power-law behaviours with different exponents. 



In this paper we have analyzed the weak lensing signatures of a fourth order (/ (R) = R n ) gravity black hole 
solution, with gravitational potential given by Eq. This introduces a new parameter a € [0, 1), which governs 
the deviation from the Newtonian gravitational potential. General relativity is contained as the special case n = 1 
(corresponding to the model parameter a = 0). A black hole in this theory generates a greater gravitational attraction 
at distances larger than r c as compared to the prediction of Newtonian gra vity. 



Lensing properties of such black holes were analyzed before in Ref. [29( , based on the small angle lens equation 
[34| . In this paper we have improved upon this approach, by employing the first order accurate Virbhadra- Ellis lens 
equation (fT6|) . 

We analyzed the dependences upon a and upon the impact parameter b of the deflection angle (|17|) . The deflection 
angle decreases with increasing impact parameter for all a. There is a transition at a critical value (b/r c ) crit = 2, below 
which the deflection angle monotonically decreases with increasing a, and above which there is a single maximum. 
This maximum value increases with the impact parameter. 

The image positions as a function of the lensing mass and source position, also the image magnifications and their 
ratio as function of source position show features similar to those in the Schwarzschild case. Nevertheless, in contrast 
with previous claims in the literature, these lensing quantities depend upon a. 

We have computed the image positions for two values of a. For the larger value of a, the image separation grows 
faster with an increase in the mass and grows more slowly as the source moves away from the optical axis. 

For the same source position the magnification factors of the images increase with <r, especially the one for the 
primary image. The increases in their ratio fii/n2 is even more significant. 

Using the most easily measurable lensing observables, the ratio of the magnifications is shown to have a power-law 
dependence on the image separations, with the power monotonically increasing with a. This behavior provides a means 
for future gravitational lensing observations to either establish the value of er or falsify the power-law gravitational 
potential discussed in this paper. Given that the next generation of radio telescopes will easily be able to resolve 
images to less than milli-arcsecond accuracy, the different rates at which the ratio of the magnifications changes should 
be able to provide a significant observational signature constraining the validity of / (i?) gravitational theories. 



IV. CONCLUDING REMARKS 
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